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We present results on the nucleon axial vector form factors Ga(^ 2 ) and G p (q 2 ) in the quenched 
theory and using two degenerate flavors of dynamical Wilson fermions for momentum transfer 
squared from about 0.1 to about 2 GeV 2 and for pion masses in the range of 380 to 600 MeV. 
We also present results on the corresponding N Xo A. axial vector transition form factors C$(q 2 ) 
and Cg(g 2 ) using, in addition to Wilson fermions, domain wall valence quarks and dynamical 
staggered sea quarks provided by the MILC collaboration. 
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1. Introduction 

Form factors are fundamental quantities which characterize important features of the hadrons 
such as their size, charge distribution and magnetization. In this work we present a calculation 
of the axial form factors of the nucleon and the nucleon (N) to A transition [p. Spontaneous 
breaking of axial symmetry in QCD is manifest in the existence of light pseudo-Goldstone bosons. 
The smallness of the pion mass is connected to the axial symmetry breaking through the partial 
conservation of axial current (PC AC) hypothesis, which relates the divergence of the axial vector 
current to the pion field. This translates into the pion-pole dominance of the induced nucleon 
pseudoscalar form factor, G p (q 2 ). Furthermore, it relates the nucleon axial charge gA, to the nNN 
coupling, gnNN via the well known Goldberger-Treiman relation (GTR), m^gA = fngnNN, where 
f n is the pion decay constant and the mass of the nucleon. Similarly, PCAC applied to the 
N to A transition, leads to the assumption of pion-pole dominance in the case of Cg and to the 
non-diagonal GTR, Im^C^ = fngnNAi where g^NA is the strong coupling constant associated with 
the 7tNA vertex. Our calculation of the nucleon axial vector form factors, G^(^ 2 ) and G p (q 2 ), as 
well as the pseudoscalar form factor G n ^, for a range of pion masses down to 380 MeV, both in 
the quenched theory and using two degenerate flavors of dynamical Wilson fermions, enables us to 
check pion-pole dominance and the Goldberger-Treiman relation. Similarly, in the case of the N to 
A(1232) weak current transition, we evaluate the corresponding dominant form factors C^(q 2 ) and 
C^(q 2 ), as well as the pseudoscalar form factor G n ^\, using besides Wilson fermions, dynamical 
staggered quark configurations generated by the MILC collaboration and domain wall fermions 
for pion masses as low as 360 MeV. The evaluation of both nucleon and N to A transition form 
factors enables us to check the q 2 -dependence of ratios of transition to nucleon form factors, such 
as G jcjva (q 2 ) / GnNN (<? 2 ) an d C5 (q 2 )/GA(q 2 ), which show a weaker quark mass dependence and are 
more likely to show less sensitivity to lattice systematics. 

2. Axial Form Factors and the Pion Pole 

We take the u- and the d- quarks to be degenerate and work in terms of the fermion isospin 
doublet y. The axial vector and pseudoscalar currents are given by A a ^(x) = y( x )7h75 \w{ x ) an d 
P a (x) = \jf(x)Y5YWi x ) respectively, where x a are Pauli-matrices acting in flavor space. The matrix 
element of the weak axial vector current between nucleon states is written as 

—u{p,s) (2.1) 

with the axial vector, Ga(<? 2 ), and the induced pseudoscalar form factor, G p (q 2 ), being functions 
of the invariant momentum transfer squared, q 2 = (p' — p) 2 . Experimentally, the axial charge 
of the nucleon, gA = Ga(0) is very well known from the neutron -decay and takes the value 
gA = 1.2695(29). The momentum dependence of Ga has been extracted from pion electropro- 
duction or quasielastic neutrino scattering experiments (see Ref. ^ for reviews) and it is con- 
ventionally parameterized by a dipole form, Ga(^ 2 ) = gA/ — fp^j > with an axial mass, Ma = 
1. 026 ± 0.0021 GeV. The induced pseudoscalar form factor G p (q 2 ) is less well studied experimen- 
tally [Q], with ordinary and radiative muon capture experiments giving different results. Both form 
factors have been studied in the context of chiral effective theories 
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In QCD the explicit breaking of axial current conservation leads to the axial Ward-Takahashi 
identity, d^A"^ = 2m q P a , where m q is the renormalized quark mass. Using PCAC we have the 
relation, d^A"^ = f n m\ii a , and therefore the pion field is expressed in terms of the pseudoscalar 
current as % a = 2m q P a / f^m^. Taking the matrix element of the pseudoscalar density between 
nucleon states, the nNN form factor is obtained via 



1 /2 

n ^xrl ' '\\i>1\m( W ( m N \ fnm\G KNN {q 2 ) , , T 3 

2m q <N(p ,s)\P \N(p,s)>= — — u{p ,s)iy 5 —u(p,s). (2.2) 

\E N {p')E N {p) I mjz-q 1 2 



Using the PCAC hypothesis together with Eq. ( |2.2[ ) we obtain the diagonal Goldberger-Treiman 
relation 

Cj tf) + ^ rC ,tf) = J- 2C "» 2 '^ . (2.3) 
4mfj 2m N m L n — q L 



In the chiral limit, the pole on the right hand size of Eq. ( J2.3| ) must be compensated by a similar 
singularity in G p {q 2 ) since G^(0) is finite. Therefore, assuming pion-pole dominance, G p {q 2 ) is 
given by 

2mfj p m\ — q 2 

and substituting in Eq. ( |2T3| ) we obtain the well known relation, m^GA{q 2 ) = fnGnNNiq 2 )- 

The invariant proton to A + weak matrix element is expressed in terms of four transition form 
factors as 

1/2 

< 



u P (p,s). (2.5) 



The form factors C^(q 2 ) and C^(q 2 ) belong to the transverse part of the axial current and, as was 
recently demonstrated [Qj, both are suppressed relative to C^{q 2 ) and C£(q 2 ). The pseudoscalar 
transition form factor, G^a. is defined similarly to G^nn, via 

2m q < A(p lS )\P \N(p,s) >= - l ^ A(pWp) ) m 2_ q 2 V(P.')^M -(2-6) 

The nNN and nNA coupling constants are defined at q 2 = nr % as guNN = GWaK^jt) anc ^ = 
G^jva The dominant form factors C^(q 2 ) and C^(q 2 ), which belong to the longitudinal part 
of the axial current, are related via PCAC to G n N&{q 2 )'- 

which is known as the non-diagonal GTR. Pion-pole dominance for C^(q 2 ), as for G p (q 2 ), leads to 
the relations C^(q 2 ) ~ m^G^A(^ 2 )A/2(m| — q 2 ) and 2rriNC^(q 2 ) = fnG n Nh{q 2 )- It is clear from 
the above relations that C5 is analogous to Ga and analogous to G„. 
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3. Lattice techniques 

The techniques used for the evaluation of the axial form factors were developed for the study of 
the electromagnetic form factors [Q, ||, |6|]. In the evaluation of the three-point functions involved in 
the calculation of the electromagnetic nucleon and N to A transition form factors we use sequential 
inversion through the sink, thereby obtaining the form factors for all momentum transfers. This 
method requires fixing the source-sink time separation, as well as the initial and final hadron states, 
but allows the insertion of any operator with arbitrary momentum at any time slice. An important 
ingredient of our method is the construction of optimal sources and sinks using a linear combination 
of interpolating fields. It was shown in Refs. ^ ||] that, for axial operators, the most symmetric 
linear combination of matrix elements that can be considered is 



£n A (0,-q ; r*;M=j) 

k=l 



c 

Am n 



(E N + m N ) (Si j + 52j + S3j) G A (Q 2 ) - 



q2 + q3)^G p (Q 2 ) 
ZmN 



(3.1) 

where j = 1,2,3 labels the spatial current direction, T projects to definite nucleon spin states and 
C is a kinematical factor. n A (p',p ;r;/j.) is the large Euclidean time limit of an appropriately con- 
structed ratio of the three-point function to two-point functions, in which time dependencies of the 
time evolution and overlaps of the initial trial state and the nucleon state cancel. Since this opti- 
mized sink, involving momentum in all spatial directions, coincides with the one calculated for the 
electromagnetic current no new sequential inversions are required for the axial vector current. The 



same holds for the pseudoscalar current and the N to A optimized sinks. As indicated in Eq. (3.1), 
we use kinematics where the final nucleon is produced at rest and therefore the momentum transfer 
q = p' — p = —p. We take — q 2 = Q 2 > with Q 2 being the Euclidean momentum transfer squared. 



# confs 


K or ami 


m % (GeV) 


M N (GeV) 


M A (GeV) 


Quenched 32 3 x 64 a" 1 = 2. 14(6) GeV 


200 
200 
200 


0.1554 
0.1558 
0.1562 
K c =0.1571 


0.563(4) 
0.490(4) 
0.411(4) 
0. 


1.267(11) 
1.190(13) 
1.109(13) 
0.938(9) 


1.470(15) 
1.425(16) 
1.382(19) 


Unquenched Wilson 24 3 x 40 [0] aT l = 2.56(10) GeV 


185 
157 


0.1575 
0.1580 


0.691(8) 
0.509(8) 


1.485(18) 
1.280(26) 


1.687(15) 
1.559(19) 


Unquenched Wilson 24 3 x 32 [§] a" 1 = 2.56(10) GeV 


200 


0.15825 
k c = 0.1585 


0.384(8) 
0. 


1.083(18) 
0.938(33) 


1.395(18) 


MILC 20 3 x 64 a" 1 = 1.58 GeV 


200 
198 
100 


0.03 
0.02 
0.01 


0.594(1) 
0.498(3) 
0.362(5) 


1.416(20) 
1.261(17) 
1.139(25) 


1.683(22) 
1.589(35) 
1.488(71) 


MILC 28 3 x 64 cT x = 1.58 GeV 


150 


0.01 


0.357(2) 


1.210(24) 


1.514(41) 



Table 1: The number of configurations, the hopping parameter, K, for the case of Wilson fermions or the 
mass of the light quarks, ;«;, for the MILC staggered quarks, and the pion, nucleon and A masses. 



4 



Nucleon and Nucleon-to-AAxial Form Factors from Lattice QCD 



Antonios Tsapalis 



In Table I we collect the parameters for our calculation. In the so called hybrid approach, 
that uses staggered sea quarks and domain wall valence quarks, we have no 6 {a) artifacts, unlike 
Wilson fermions where cutoff effects are G(a). We use the same parameters for the domain wall 
operator as those used in Ref. [§], namely we take the length of the fifth dimension, L^ja = 16 and 
the valence quark mass that is tuned to reproduce the pion mass calculated with the staggered quark 
action. Finite volume effects can be assessed by comparing results at the lowest pion mass. It was 
shown [|l|] that results on the 20 3 spatial volume were consistent with the results on the 28 3 lattice. 
The source-sink separation is optimized so that, on the one hand, it is sufficiently large to ensure 
ground state dominance and, on the other, small enough so that gauge noise is kept at a minimum. 
In all cases, we use the non-perturbatively determined value for the axial renormalization constant, 
Z A . Note that the pseudoscalar renormalization constant is not needed for the quantities under 
consideration here. 



4. Results 



We first consider the ratios G nNh {Q 2 )/G KNN {Q 2 ) and 2C^{Q 2 )/G A {Q 2 ) 



We note that, in 

the ratio G JZ N^(Q 2 )/G n ^^(Q 2 ), the renormalized quark mass cancels eliminating one source of 
systematic error. 
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Figure 1: The 



GWa (Q 2 ) /G n NN {Q 2 ) 



ratio of form factors 
as a function of Q 2 
for Wilson fermions for the quenched theory, 
denoted by Nf = 0, and for two dynamical Wilson 
quarks, denoted by Nf = 2. 



Figure 2: The ratio of 2C%(Q 2 )/G A {Q 2 ) as 
tion of Q 2 . The notation is the same as 
Fig. 1. 



a func- 
that of 



As can be seen in Figs. [I] and ||, these ratios show weak dependence on the quark mass and 
are therefore more suited for comparison with physical results. Both these ratios are Q 2 inde- 
pendent and can be thus fitted to a constant. Fitting the quenched data, which carry the small- 
est statistical errors, we obtain the value of 1.60(2) shown by the dashed line. Taking the ratio 
of the diagonal and non-diagonal GTRs m^Ga = /kG^nn and 2m^C^ = fnGnNk, we find that 
G nN ^Q 2 )lG nNN {Q 2 ) = 2C*(Q 2 )/G A (Q 2 ). In Fig. |, we show the ratio 2Cj{Q 2 ) /G A (Q 2 ), which 
is indeed also Q 2 independent, and fitting to the quenched data we find the value of 1.63(1) shown 
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by the dashed line. They are also approximately equal in the unquenched case [|I|]. Therefore, on 
the level of ratios, the GTRs are satisfied. 
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Figure 3: G A {Q 2 ) (upper) and G p (g 2 )(lower) for 
the quenched theory and for two dynamical Wilson 
fermions. Results in the hybrid approach are from M\. 



Figure 4: C^(Q 2 ) (upper) and C^(Q 2 ) (lower) 
for Wilson fermions and in the hybrid approach. 



In Fig. ^J, we show our results for the nucleon axial form factors using Wilson fermions. We 
also include recent results obtained in the hybrid approach [[J using similar parameters to those 
used in our evaluation of the N to A axial form factors. The main observation is that, at the smallest 
domain wall quark mass, results in the hybrid approach deviate from quenched results. In partic- 
ular, we note that the value of the nucleon axial charge gA becomes larger in the hybrid scheme 
approaching the experimental value. The dotted line shows the dipole fit to the experimental data. 
As can be seen, lattice results fall off slower than experiment. However, the large deviations ob- 
served at low Q 2 for a pion mass of about 360 MeV point to large pion cloud effects, which tend to 
increase the form factors at small Q 2 . The solid line is a dipole fit of the quenched Ga(Q 2 ) at a pion 
mass of 41 1 MeV, yielding an axial mass Ma = 1 .58(3) GeV. We note that a fit to an exponential of 
the form goexp(— G 2 a l so describes the lattice data, yielding a curve that is indistinguishable 
from the dipole fit. 

Very similar behavior is observed for the transition form factor Cj (Q 2 ) shown in Fig. [l], where 
a dipole fit, shown with the solid line, provides a good description of the lattice results. Again the 
axial mass obtained by fitting the quenched results at the lowest mass is larger than the value of 



1.28(10) GeV extracted from available experimental results [10]. A fit to an exponential form also 
provides a good fit to the lattice data. 

Having fitted G A {Q 2 ) and C*(Q 2 ), the £ 2 -dependence for the form factors G P (Q 2 ) and C^(Q 2 ) 
can be obtained using pion-pole dominance. The resulting curves are shown by the dashed lines in 
Figs. U and |] for the quenched theory at the lowest pion mass and show deviations at low Q 2 . The 
dotted line is the corresponding result in the hybrid approach at the smallest pion mass. In addition 
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we show curves that are obtained using 

G P m _ g o (41) 



G A (Q 2 ) (Q 2 /m 2 + \) 

3< 



with a corresponding expression for C^(Q 2 )/C^(Q 2 ), where go and m are treated as fit parameters. 



As expected, this provides a good description of the <2 2 -dependence for G P (Q 2 ) and Cg (Q 2 ) as 
shown by the solid lines, which correspond to the fits of the quenched data at the lowest pion mass. 

5. Conclusions 

We present results for the nucleon axial vector form factors G A (Q 2 ) and G p (Q 2 ), as well as 
for the corresponding N to A axial transition form factors C^(Q 2 ) and Cg (<2 2 )- The nNN and nNA 
form factors G n NN(Q 2 ) and G^jva(8 2 ) are also evaluated. One of the main conclusions is that G n NN 
and Gjzna have the same Q 2 dependence yielding a ratio of G^a(2 2 )/G^ww(2 2 ) = 1.60(2) in 
good agreement with what is expected phenomenologically. The ratio 2C^(Q 2 )/Ga{Q 2 ) = 1.63(1) 
is also independent of Q 2 . Equality of these two ratios implies the Goldberger-Treiman relations. 
We also studied the <2 2 -dependence of the form factors separately. Comparing quenched and un- 
quenched results at pion mass of about 360 MeV, we observe large unquenching effects on the low 
g 2 -dependence of the four form factors, G A (Q 2 ), G P (Q 2 ), Cj{Q 2 ) and C%(Q 2 ). This confirms the 
expectation that pion cloud effects are expected to be large at low Q 2 . Further study of pion cloud 
effects using lighter quark masses on a finer lattice is called for. 
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